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Numerical relativity, applied to collisions of black holes, starts with initial data for black holes 
already in each other's strong field. For the initial data to be astrophysically meaningful, it must 
approximately represent conditions that evolved from holes originally at large separation. The initial 
hypersurface data typically used for computation is based on mathematically simplifying prescrip- 
tions, such as conformal flatness of the 3-geometry and longitudinality of the extrinsic curvature. 
In the case of head on collisions of equal mass holes, there is evidence that such prescriptions work 
reasonably well, but it is not clear why, or whether this success is more generally valid. Here we 
study these questions by considering the "particle limit" for head on collisions of nonspinning holes, 
i.e., the limit of an extreme ratio of hole masses. The mass of the small hole is considered to be a 
' perturbation of the Schwarzschild spacetime of the larger hole, Einstein's equations are linearized 

0^ . in this perturbation, and described by a single gauge invariant spacetime function ij), for each mul- 

0^ ' tipole. The resulting quadrupole equations have been solved by numerical evolution for collisions 

starting from various initial separations, and the evolution is studied on a sequence of hypersurfaces. 
In particular, we extract hypersurface data, that is t/j and its time derivative, on surfaces of constant 
' background Schwarzschild time. These evolved data can then be compared with "prescribed" data, 

evolved data can be replaced by prescribed data on any hypersurface, and evolved further forward 
in time, a gauge invariant measure of deviation from conformal flatness can be evaluated, and other 
' comparisons can be made. The main findings of this study are: (i) For holes of unequal mass the 

use of prescribed data on late hypersurfaces is not successful, (ii) The failure is likely due to the 
inability of the prescribed data to represent the near field of the smaller hole, (iii) The discrepancy 
in the extrinsic curvature is more important than in the 3-geometry. (iv) The use of the more general 
conformally fiat longitudinal data does not notably improve this picture. 
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I. INTRODUCTION AND BACKGROUND 



Q^. The collision of black holes is of great interest as both astrophysics and as a strong-field gravitational interaction with 
^ ' no Newtonian analog. Black hole collisions also may provide the strongest, most observable, source of gravitational 
radiation that can be detected by gravitational wave observatories now under construction |0] . The strong field nature 
■ of the process means that the nonlinear character of Einstein's theory plays a crucial role, so that estimates based on 
. !^ \ linearized theory are unreliable. For the last stage of coalescence of two rotating holes in a decaying orbit, we have 
only dimensional estimates. The timeliness of the problem has given rise to an effort to attack the problem by solving 
5h Einstein's equations numerically on supercomputers Q| . 

^ The numerical solution starts with some "initial value" data, a solution of a subset of Einstein's equations on an 

initial hypersurface, and then evolves this solution forward in time. The specification of the initial data is tantamount 
to specifying what problem it is that the computer will be solving, so the initial data must encode the physical 
parameters (mass, spin, location, momentum) of the colliding objects. These physical parameters are meaningful 
and well understood for isolated holes, but become increasingly ambiguous as the initial separation between the 
holes decreases, and the field of each hole has a strong effect on the spacetime geometry of the other. Ultimately, at 
sufficiently small initial separation, it becomes impossible to make any clear statements about the physical parameters 
of an individual hole. 

In principle, one would like to avoid this difficulty by specifying initial data on a hypersurface in the distant past 
when the holes were separated by a distance many times the horizon size. When the initial infiuence of each hole 
on the other becomes small, the uncertainties in the choice of initial data become negligible. Unfortunately, it is 
not possible to start the numerical evolution a long time before the collision. Einstein's equations are a nonlinear 
hyperbolic system, and for such systems instabilities in numerical evolution seem to be a generic feature. While it 
can be expected that some progress will be made in improving the codes and suppressing the instabilities, it is not 
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likely that in the near future these codes will be able to evolve initial data for more than a few dynamical times, the 
characteristic time scale for black hole processes (around 10~^sec for a solar mass hole). This limitation precludes 
starting the dynamics with the holes at large distances and means that the initial data must be given on a hypersurface 
on which the hole-hole interactions are already strong. If the resulting numerical evolutions are to be astrophysically 
meaningful, it is crucial to have a way of connecting data on this "late" hypersurface to an astrophysical precursor 
configuration with the holes interacting weakly. In other words, to understand the meaning of initial data on a 
"late" hypersurface, we must know where that data came from. One approach to providing this connection is to 
use approximations such as higher-order post-Newtonian methods in order to understand the evolution of the system 
from an early stage to a stage at which the holes are interacting with intermediate strength. Such an approach may 
turn out to be sufficient if the numerical evolutions can be stabilized for a relatively long time. 

There is reason to hope that it might be possible to give initial data at late times, with relatively little difficulty. The 
reason is based on the success of a "two phase" approximation method used by Abrahams and Cook ^ and further 
investigated by Baker and Li These studies dealt with the head on (zero angular momentum) collisions of two 
equal-mass nonrotating holes. Data were specified on a very late hypersurface, simply by using a "standard" initial 
value solution formally representing equal mass holes moving towards each other. This standard solution contained 
a parameter P which, in the case of an isolated hole, agreed with the momentum of the hole. On the very late 
hypersurface the separation and momentum were set at the values dictated by Newtonian gravity theory if the holes 
had started infall from some large initial separation. The standard initial value data was evolved forward in time (with 
the use of an approximation method). The radiated energy found by this approximation could be compared with 
the "correct" radiated energy, since numerical relativity results were available for such collisions. The approximation 
method agreed with numerical relativity to rather good accuracy. This is somewhat a surprise since the "standard" 
numerical relativity initial data should be very different from the "correct" evolved data. The implication of this 
agreement is that results for gravitational radiation might not be sensitive to all details of initial data. 

We explore this question in the present paper by comparing prescriptions for imposing data on hypersurfaces with 
the "correct" data evolved from an earlier configuration. In order to do this, of course, we need the ability to evolve 
a solution for a relatively long time, something that cannot yet be done with the nonlinear Einstein equations. We 
choose therefore to use the particle approximation, the approximation in which the mass m one of our holes, is much 
smaller than the mass of the other M . We then treat the low- mass hole as a particle of negligible size, and we 
treat m/M as an expansion parameter. By doing first-order perturbation theory with this parameter, our evolution 
equations become linear, so that a stable numerical scheme is straightforward to develop. 

In an earlier paper |^ (Paper I), we developed the basic mathematics of the particle limit approximation, and of 
evolving data that was initially stationary. In that paper, Laplace transform methods were used. The method for the 
present results is a direct numerical solution of the partial differential equations in radius and time (after multipole 
decomposition). We have compared results of the two completely different methods, and found agreement within the 
estimated numerical accuracy, thereby adding confidence that there are no mistakes in our numerical methods. 

The paper is organized as follows. In Sec. II we give the mathematical formulation of the problem. We base this 
formulation on the gauge- invariant Moncrief |^ approach to perturbations, so that the variable we use is invariant 
with respect to perturbative changes in slicing. We need only specify the slicing to zero order in the perturbation, and 
to that order we take our hypersurfaces to be slices of constant Schwarzschild time. In Sec. II we discuss in detail the 
freedom in initial data and the standard choices made for reducing that freedom. We go on to explain the particular 
physically-motivated choice that we make for arriving at a definitive prescription of initial data. In Sec. Ill, we start 
by explaining our numerical method for evolving ij). The scheme for solving the linear hyperbolic equation obeyed by 
ij) is straightforward. The only issue that deserves attention is the source term in the equation for ip. In the particle 
limit, the small hole is treated as a point particle, and this particle is associated with a stress energy source that 
enters the tp equation as the derivative of a delta function of radial position. We describe in Sec. Ill how the singular 
source is handled. Numerical results are presented and compared in Sec. III.B. In Sec. IV we discuss the results and 
their meaning. 

II. MATHEMATICAL FORMULATION 

A. Moncrief-Zerilli formalism 

We will describe the perturbations to the Schwarzschild background, due to the particle, in the notation of Regge 
and Wheeler Q (hereafter RW). The symmetry of the straightline infall of the particle means that there will be no 
odd parity perturbations. For the even parity perturbations the general form of the line element, with perturbations 
of a specific multipole index is 



2 



ds^ = dsl + (1 - 2M/r){H^Ym)dt'^ + (1 - 2M/r)-\H^Yio)dr'^ 
+r^{K'^Yeo + G^d^Ym/de^)de^ + r^{sm^ OK'^Ym + sinO cosedYio/d0)d^^ 
+2H[Ymdtdr + 2h^o{dYeo/de)dtde + 2h[{dYm/de)drde . 



(2.1) 



Here dsg the unperturbed line element for a Schwarzschild spacetime of mass M, where H^, Hf, H2, hf), h{, and 
are functions of r, t, and Ifo(^) are the m — spherical harmonics. For simplicity, from here on we shall drop the 
£ index on perturbation functions. 

There are two closely related formalisms available for describing the evolution of the perturbations in terms of a 
single wave function and a single wave equation. The method due to Moncrief [D uses only information {H2,hi,K, G) 
about the 3-geometry of a i= constant hypersurface, and is gauge invariant, i.e., independent of the choice of lapse 
and shift, and of diffeomorphisms on the hypersurface. The second method, due to Zerilli |8|, relies on a specific 
gauge, the gauge choice introduced by Regge and Wheeler in which hQ,hi, and G are set to zero. 

As in Paper I we shall use the Moncrief formalism. The Moncrief wave function ip, in terms of the perturbations in 
the RW notation, is 



V'(r,t) = 



A + 1 



K 



r - 2M 



{H2 - rdK/dr} 



Xr + 3M 

where we have used Zerilli's normalization for tp and his notation 

A= (£ + 2)(^-l)/2 



+ 2(1 - 2M/r) {r^dG/dr ~ 2hi) 



(2.2) 



(2.3) 



For simplicity, we shall occasionally present some equations restricted to the RW gauge. Since the Moncrief wave 
function is gauge invariant, this gauge restricti on a llows simplicity of presentation, with no other consequences. In 
this gauge, we omit the G and hi terms in Eq. (2.2). 

The basic wave equation for an infalling particle is given in Paper I as 



(2.4) 



Here r* = r + 2AI \n{r /2M — 1) is the Regge- Wheeler [[7[ "tortoise" coordinate and Ve is the Zerilli potential (given, 
e.g., in Paper I). For a point particle of proper mass mg, the stress energy is given by 



(2.5) 



where is the particle 4- velocity. The two dimensional delta function 5^ {Q) gives the angular location of the particle 
trajectory 



5\^)=Y,Y,,r,{e,cP)Yl, 



^r,o(0)V(2^+l)/47r 



(2.6) 



with the last expression applying for infall along the positive z axis. The time dependent location of the particle rp (t) 
follows from the geodesic equation and, for a particle starting from rest at t — 0,r = tq, is the inverse of 



2M 



ro \2MJ \2MI 



1/2 



1 



AM / ro \^/^ / 2M 

~'\2m) 



3/2 



2M 



■ arctan 



^-1 



In 



1 - (1 - AM/ro) {rp/2M) + 2^1 - 2M/ro^rpl2M^\ ~ rp/rp 

(rp/2M) - 1 



(2.7) 



From the particle stress energy in (2.5), the source term Si{r^t) on the right hand side of (2.4), is given in Paper I 



as 



Si{r,t) 



2(1 - 2M/r)K 



r\l-2Mlr)—5'{r-rp[ 



r(A + l)(Ar + 3Af) 
r(A + 1) - M 3M[/"r(l - 2M/rf 



2[/o 



Ar + 3Af 



5{T-rp[t]) 



(2., 
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where 



K EE 87rmov/(2£+ l)/47r . 
The total radiated energy after a given to is 

1 (^ + 2)! 



Energy = 



647r - 2) 



dt 



(2.9) 



(2.10) 



B. Hypersurface data for t/j 



We choose our data, on a hypersurface of constant t, to correspond to the choice usually made for work in numerical 
relativity. The data, for both the initial 3-geometry and extrinsic curvature, are based on the work of York and 
collaborators. This approach takes the 3-geometry of the initial hypersurface to be given by da^ = ^'^da\^^. The 
conformal factor $ must, in vacuum, satisfy the equation 
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(2.11) 



where is the Laplacian taken with respect to the flat background, and where Kij is the extrinsic curvature. 
For a constant t slice of the unperturbed Schwarzschild spacetime, the extrinsic curvature vanishes. It follows that 



Kij is perturbative, and hence the right hand side of (2.11) is a higher order perturbation, which, in our first order 



perturbation calculation, we can ignore. We denote coordinates in the conformally related flat space with bars, such 
as r, z, and we note that $ must be a harmonic function of these coordinates. We make two different choices for $. 
The flrst is the solution discussed by Brill and Lindquist |l^, and which we will call a BL type solution. This solution 
follows if we treat 2$ as if it were the Newtonian potential of two points, one of mass M, and one of mass mp. If the 
first mass is located at the origin of the flat coordinates, and the second is at (flat) coordinate position z = Zp, then 
the BL solution takes the form 



$BL = 1 



M 
2f 



mp/2 



(2.12) 



•ipcz I 



ZpGz are coordinate singularities, they are not geometric singularities. Near 



Though the origin and the point fe. 
these coordinate points the divergence of the conformal factor means that the "points" are, in fact, asymptotically 
flat regions, so the geometry given by ( ^.12 ) is actually that of three asymptotically flat spaces connected by two 
throats. For our purposes here, we take the "particle limit," the limit that rup <^ M . 

The BL solution is the only choice that was considered in Paper I. A second initial conformal factor that can be 
considered is that which is "reflection symmetric," in the sense that there are only two asymptotically flat regions, and 
they are isometric to each other. The method of constructing such solutions was given by Misner ||l^, and requires 
that "image" points be used in the conformal space, in such a way that the geometry is symmetric with respect to 
inversion through a sphere about each of the point singularities. Here we make the solution symmetric only for a 
sphere about the origin. The singularity at fe^ = Zpez, in the limit, will be taken to have a zero size inversion sphere 
around it (the particle limit) and the images inside it will be meaningless. (Alternatively one can construct a solution 
reflection symmetric for two spheres, and subsequently take the limit; the result is the same.) The conformal factor, 
for this "Misner" case, is 



1 M 

2r re 



TOp/2 



gc/2 



(2.13) 



It is straightforward to verify that this Misner solution is symmetric for inversion about r — AI/2 if the parameters 
are chosen to be milage ~ rUpM / {2zp) and Zimagc — {M /2)'^ / Zp. 



^imagc, •^linage 



For either the BL form or the Misner form the line element for the 3-geometry can now be written 



da' 



l + M/2f+ ai{f)PeXcos9) 

1=0,1,2, ■■■ 



{df^ + f^dif) 



(2.14) 
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with dQ"^ = dO"^ + sin^ Bdcj)^ . The requirement that $ be harmonic means that must have the form a^f^ + bef~^^~^^'> . 
We must now put the metric into a form for comparison with the Schwarzschild metric. To do this we introduce a 
Schwarzschild-Uke radial coordinate r related to the flat space coordinate f by 



(2.15) 



This is only one of many possibilities for such a transformation, but it is the si mplest for computation, and has been 
very convenient for perturbation analysis This transformation puts (2.14) into the form 



da^ = 



1 + 



2mp/r 
+ M/{2f) 



J2 Mf)p,{cose) 



«=0,1,2--- 



dr2 



1 - 2M/r 



^dn^ 



(2.16) 



where we have kept only terms first order in m^, and, as in Paper I, have introduced the notation ^-^(r) = {2r / mp)ai{r) 
The perturbations, in the RW notation, can now be written 



2mp/r 
1 + M/2f 



2e + i 



(2.17) 



with hi = G = 0. These perturbations of the initial hypersurface turn out to be in the RW gauge, so we can compare 
them with the hamiltonian constraint, in the RW gauge, as given by Zerilli p3[ |, and in Paper I, for a particle of 
proper mass uiq. 



1 



2M\ d^K 
r 



5M\ IdK 



1 



2M\ ldH2 
r ) r dr 



dr 



-87rmr 



2r2 



{K + H2 



1 



[H2 - K) 



21- 



(l - ls(r 



T r' 



5{r~rp) , 



(2.18) 



in which the right hand side comes directly from the stress energy expression (2.5). 

For the choice of a conformally flat 3-geometry we have H2 = K and Eq. ( 2.18]) becomes a differential equation for 
K. A particular solution of this equation, well behaved at infinity, is 



j^partic 




r > Zp 

f < Zr, 



(2.19) 



This solution, in fact, corresponds precisely the BL solution of Eq. (2.12). For Zimagc < M/2 and r > M/2, a 
homogeneous solution well behaved at infinity is 



^homog 



Air 



2m\, 



21+1 V 1 + M/2f 



^ 2^imagc ^ 



(2.20) 



If we take z\ 



image 



{M/2Y/zp and m 



image 



solution due to an "image" mass, inside the horizon sA r = z 



nipM/ {2zp), then the homogeneous solution can be considered to be the 



image ; 



and it is straightforward to verify that with these 
choices of ^image and mimage the Misner solution is K'^^^ = if P^ftic _|_ ^homog^ There are, of course, other possibilities. 
We could, for instance, consider, with the same choices of Zimage and milage, the combination K""^^' = i5:partic_;;^homog^ 
For this solution K and "0 would vanish on the horizon. This property turns out to be preserved by evolution as is 
made clear in the discussion below. 

By using (2. IE) and integrating across the singularity in (2.18) at r = = Zp[l + M/(2zp)]^, wc find that 



AK,r = dK/drl^^+ - dK/drl^^ 



-8T:m 



^{2 e+i)/ATT 

'P-1/2 3/2 



Zp'r;''y/l~2M/rp 



(2.21) 



The 4-velocity component JJ" is given by 
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so it follows from ( 2.21 ) that 



U 



_ v/1 - 2M/ro 



1 - 2M/rp 



(2.22) 



'p 



1 / / 2A/ 

= ^mo 1 + Wl 



' 1 - 2M/ro 
1 - 2M/rp 



(2.23) 



Note that the relationship of the mass parameter rrip and the particle mass mo is the same for the Misner case as 
for the BL case, and hence the same as in Paper I. This must be true, of course, since the symmetrizing image term 
addition in ( ^.20 ) has no discontinuity at r = rp. 

The initial hypersurface, on which the particle and the hypersurface data are stationary, is always denoted by 
t = 0. On hypersurfaces with t > we shall limit our choice of prescribed data to conformally flat 3-geometries. The 
choice of conformal flatness is not, however, preserved by evolution. This means that the 3-geometry on constant t 
hypersurfaces will in general not be conformally flat. On such a hypersurface, we will have numerical values only for 
a single function, the Moncrief function tp. It turns out that we can test ip for underlying conformal flatness. The 
combination of metric perturbations 



(2.24) 



is gauge invariant, and clearly vanishes for a 3-geometry that is in conformally flat form, with hi = 0, G = and 
H2 ~ K. [See Eq. (2.1).] The computation of this gauge invariant quantity from ip{r) is most easily described in the 
RW gauge, where /conf reduces to H2 ~ K. From Eq. ( p. 2D in the RW gauge, and from Eq. ( 2.18| ), which is already 
in the RW gauge, it follows that (for r ^r^) 



6M2 + 3MAr + A(A + l)r2 , 
= ^TT h — w 



1 - 



2M\ dip 
dr 



From this result for K, and from Eq. (2_^) in the RW gauge, one flnds H2 to be (for r ^ 



' " r2(Ar + 3Af)2 



3AP - XMr + Xr 
r{Xr + 3M) 



-drip +{r- 2M)dliP 



(2.25) 



(2.26) 



The difference between the expressions in Eqs. (IIH) and ( |2^ ) gives /conf, the gauge invariant measure of the 
deviation of the 3-geometry from conformal flatness. 



C. Hypersurface data for ip 



The integration of the wave equation (2.4), to evolve forward in time from a, t = constant hypersurface, requires 
that we specify both ip and ip at the hypersurface. In Paper I our initial hypersurface was one of time symmetry, so 
we had ^ = initially. We now take the particles to be moving on the hypersurface, and for the hypersurface value of 
the extrinsic curvature we follow the prescription of Bowen and York [ll4| . In that prescription one chooses a slicing 
such that the trace of the extrinsic curvature vanishes, and one deflnes a quantity Kij, related to the extrinsic 
curvature Kij by 

= , (2.27) 

where 7^^^* is the flat 3-metric. The curvature measure Kij is considered to be a tensor in the conformally related 
flat space; its indices are raised and lowered with ^ff'^- The momentum constraint turns out to be equivalent to the 
requirement that Viif*-' vanish outside the sources, which we take to be points in the conformally flat 3-geometry. 
The Bowen- York prescription for initial data is that Kij be "longitudinal," i.e., derivable from a vector Wi according 
to 

K,, = V,W, + V, Wi - (2/3)7«fVfcW'= . (2.28) 
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Our choice to describe a moving hole is a solution closely related to those used in most numerical relativity 
simulations. To describe these solutions we use the geometry of the conformally related flat space pictured in Fig. 
The symmetry axis is taken to be the z coordinate axis and the perturbative particle is located at z = Zp. The vectors 
f and p are the displacements (in the flat space) to a field point, respectively from the coordinate origin and from 
the location of the moving "particle." The unit vector n, defined to be p/|p| is the direction to the field point from 
the moving particle. In terms of the notation in this figure, we define two solutions of the momentum constraint 
S/iK^^ — 0. The first, which we shall call the basic solution, /.fbasic^ 



basic 



2p^ 



This is a longitudinal solution that corresponds to the vector 



= -jz[7P^+n,P''nk] . 



(2.29) 



(2.30) 



The parameter P is the momentum of a moving hole in the conformally related flat space. For our particle, P is 
related to the true 3-momentum (i.e., the spatial part of the 4- mome ntum) p by \P\ — ^^\p\. Here $ is the conformal 
factor and, to lowest order, has the form 1 + M/2f, as in Eq. (2.16). In our head on collision the only nonvanishing 
component of momentum is the radial component, so we have 



P=\P\^ $V/Vl-2M/r, 



(2.31) 



with = moU^ the usual contravariant radial component of the 4-momentum in the Sc hwar zschild background 
coordinates. This relationship can be verified by using the extrinsic curvature given by Eq. ( 2.29| ) in the momentum 
constraint Rtr = SirTtr and taki ng th e right hand side to be the stress energy of a radially moving particle. 

The "basic" solution in Eq. ( |2.29| ) is the simplest solution of the momentum constraints. It is straightforward, 
though tedious, to decompose the conformal extrinsic curvature into spherical harmonics. To express the results we 
adopt the following notation, analogous to that of Regge and Wheeler [Q for the metric perturbations: 



Kr 

Kf. 



^ Ai{r;Zp)Yio 
= 1,2, ■■■ 

Be{f;zp)dY,o/de 

= 1,2, ■■■ 

2 J2 [lCdf;zp)Yeo + ge{f;zp)dYio/d''0^] 

£=1,2, ■■■ 

Hin^e J2 [ICi{f;zp)Ym + cote ge{f-Zp)dYeo/de] 



(2.32) 
(2.33) 
(2.34) 
(2.35) 



Other components vanish due to the azimuthal symmetry. 

For the expansion of the basic extrinsic curvature, the explicit expressions for the coefficient functions are 



y^basic 
^basic 
/^basic 
i^basic 



An 



2£+l \ 2f2 



P 



An 

2JT1 V 2f 

47r / P 

2e + i V2f2 

47r / P 



J-i 



{l + l){l + 2) 



I 



21- 1 



jl - 6) 
2/ + 3 



{1-2) ^,{l-6) 



21- 1 

51 
21 - 1 



q 



21 + 3 
2(^-6) 



2/ + 3 



2£+l \ 2f2 



J-i 



(1-8) , 2(^-6) 



21 - 1 



21 + 3 



(2.36) 
(2.37) 
(2.38) 
(2.39) 



Here q = Zp/r < 1 and for q > 1, the above expressions are valid if £ is replaced throughout by —£ — 1 everywhere on 
the righthand side, except in the normalization constant ^/47ij2r+T. 

The parameter P is taken to represent the (perturbative) momentum of the particle. We take our asymptotically 
flat coordinate system to be one in which there is no net momentum. This means that the large "background" hole 
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must have momentum which is in some sense equal in magnitude to the momentum of the particle. If we take the 
origin of our coordinate system, in the conformally flat space, to be (in some sense) at the center of mass, then the 
coordinate singularity representing the background hole is at a coordinate distance Zbh — — {mp/M)zp from the origin. 
We can get the contribution to the extrinsic curvature from the background hole by reversing the sign of P in the 
expressions above, and by substitu ting Zb /i in place of Zp. But Zbh is a first order quantity, so (zwi)" terms can be 
kept only when n = 0. From Eq. ( 2.36 )-( 2.39 ) we see that this means that only £ = 1 terms can be kept. But the 
£ = 1 even parity terms have no physical content, and are not coupled to gravitational radiation. In treating the 
perturbative extrinsic curvature, therefore, we consider only the contributions from the particle. 

To use the above results in calculations, we must relate the extrinsic curvature to the Moncrief wave function ip. It 
is useful to have at hand the relations between the extrinsic curvature and the metric perturbations. This is found by 
considering at— constant slice of the perturbed spacetime. It is straightforward to calculate the extrinsic curvature 
of this slice in terms of the perturbations of the spacetime metric. In the RW notation (but not the RW gauge), the 
results are 

(2.40) 
(2.41) 
(2.42) 
(2.43) 



A 


= $ 


-2 


(1 - 2M/r)-^/^H[ - 


^(1- 


^ 2Mlr)-^/^H2 ^ 


-(l-2M/r)-3/2(A//r2)iJi 


B 


= $ 


-2i(l-2M/r)-i/2 


Hi 


+ h'o- 


- hi — 2ho/r 






JC 


= ^' 


-2 


r(l - 2M/rfl'^ 


Hi - 




- 2M/r)-^^^k 




Q 


= $ 


-^(l-2Af/r)-i/2 


ho - 




1 







where prime indicates derivate with respect to r, and dot with respect to t. 

To find the initial value of ■(/', we use the prescription given by Abrahams and Price |lE 



and we write 



(2.44) 



in which "Q" represents Moncrief combination ||^ of perturbations of the metric Sgij, and of the r derivative of Sgij. 
The initial value of ■0 is then the initial value of 

i^^-2Qi{l-2M/r)^/^6K^,,(^{l~2M/r)^/^SK^jy\ . (2.45) 



By comparing Eqs. (|2T|) and (|2.32|)-(|2.35|) we find that the substitutions needed are: 

(1 - 2A-//r)"^' 



'H2 

K 
G 
hi 



{df/drf v/l - 2M/r 
-^K. {f/rf - 2M/r 



[f/ry yJ\^^2Mlr 
$"^6 (f/r) - 2Af/r. 



(2.46) 
(2.47) 
(2.48) 
(2.49) 



With these substitutions made in Eq. (2.2) we get a specific expression for ^, as a function of r, on the hypersurface 

with the express 

\h2- rdkjdr^ 



An independent derivation of the result starts with the expression 

- 2M 



A + 1 



K 



Xr + 3M 



+ (r - 2M) r^dG/dr - 2h_ 



(2.50) 



where dots indicate d/dt. With Eqs. (2.40)-(2.43). the right hand side of Eq. ( 2.5C ) ca n be reexpressed in terms of 
A,B,JC,G- We have checked that the result is identical to that found from Eqs. (2.45)-(2.4£), i.e. 



Tpir, t) 



2f^\ - 2M/r . „2 



(A + 1)7-2 



^/l - 2M/r 
(2A + 6M/r) 



-4(A+l)(f2g) + 



M 



+ 



(r-M) 



(2.51) 



In describing the general choice of '0 on a constant t hypersurface, it is useful to review the equivalent question for 
ifj. If we restrict our choice to a conformally flat 3-geometry, then the hamiltonian constraint, and the constraint of 
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conformal flatn ess, r educes the specification of the initial geometry, for each multipole, to a single ordinary differential 
equation, Eq. ( 2.18|) with H2 = K. The on ly ch oice to be made is the constant, specifying how much of the 
physically well-behaved homogeneous solution (2.20) is to be added to a particular solution; specifyin g this constant 
is equivalent to choosing between, for exam ple, t he nonsymmetric (Brill-Lindquist) solution of Eq. ( 2.12 ) and the 
reflection symmetric (Misner) solution of Eq. (2.13). If, on the other hand, we do not restrict ourselves to a conformally 
flat 3-geometry, we have a functional degree of freedom in the initial geometry. (Equivalently, in the RW gauge we 
have the freedom to specify H2 — K for each mul tipole.) Whether or not we restrict to a conformally flat 3-geometry, 
the Moncrief wavefunction -0, deflned in Eq. (2.2), is fixed once the 3-geometry is specified. Furthermore, ■(/;, is totally 
gauge invariant. 

For -0 the situation is closely parallel. Analogous to the condition on 3-geonietry, that it be conformally fiat, we 
now have the special condition that the extrinsic curvature be purely "longitudinal," i.e., derivable from a potential 



according to Eqs . ( 2.27 ) and ( 2.2g ) . For our axisymmetric infall, the mathematical problem of finding the Wi potential, 
for each multipole, is the problem of finding the two components Wp and Wg. These can be expanded in multipoles 
and written as: 



(2.52) 



The nontrivial momentum constraints, equivalent to the t9 and tr components of the einstein equations, then constitute 
two second order differential equations for Ai{f) and Bi{f). A particular solution to the momentum constraints is, of 
course, given by the fundamental solution in Eq. ( |2.29D . This solution corresponds to the multipoles 



^basic / 



47r P 



2i+l Af 



21+3 ^ 



(2.53) 



= (r-) 



47r P 



2£-(-l 4 



(8-^) .-1 , (^-6) 



2^-1 



2^ + 3 



q 



(2.54) 



The above expressions are valid only for q = Zp/f < 1. For q > 1 the replacement i — > —£ — 1 must be made inside 
the square brackets. 

One finds that there are only two solutions well behaved at infinity, and that all solutions are well behaved at the 
horizon. There are then two we ll-beh aved homo geneous solutions of the momentum constraints that can be added to 
the particular solution of Eqs. ( 2.53 ) and ( 2.54 ). 

Homogeneous longitudinal solution to the momentum constraints can be found directly from the differential equa- 
tions. Alternatively one can con struct those from known solutions. One such solution is an "image solution," the 
basic solution given in Eq. (2.30), but with Zp replaced by the location ^imago < M/2, so that the position of the 
z = ^imago singularity is inside the horizon of the background hole. (If the choice Zjmagc — iM/2) /zp is made, then 
the singularity is the "image" of the particle under inversi on th rough the spherical surface at f = M/2.) The image 
solution is, then, the following simple modification of Eq. (2.53) and (2.54): 



1 image /-N _ 



(f) 



An P 



21+1 if 



e{£ + 7) , {£+i)(6-i) 



21 ~ I "^'"lago 



,t+l 



2£ + 3 



(2.55) 



B 



image /-N _ 



(f) 



47r P 



21+1 4 



(8-^) £-1 (^-6) ,+1 

2£_l '^^imago '^21+3 



(2.56) 



Here q. 



image 



^image 



/r and is always less than unity for r outside the horizon (i.e., r > M/2). Similarly, we ca n find 



t he mu ltipole decomposition of the corresponding extrinsic curvature, simply by replacing q by Qimage in Eqs. ( 2.36| ) 

-(EID- __ 

A second longitudinal solution can be constructed from the "a" solution given by Bowen and York [Ijjl^] and 
elsewhere: 



if" 



— [P,n, + P,n, + (7ff - hn^n,) P'^n^] 



(2.57) 



The equivalent Wi vector is 
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(2.58) 



This solution has been of interest in connection with inversion symmetry , but we us e it h ere with no prejudice about 
symmetry of the solution under inversion. To make the solution of Eqs. (2.57) and ( ^.58| ) a homogeneous solution we 
need only place the p = singularity inside the horizon, say on the z axis at Zg < M/2 . 
The multipole decomposition of the a solution, in the notation of Eqs. ( 2.321) - ( 2.35 ) is 



A? 



An P 



21 + 1 2f4 



An P 



21 + 1 2f 3 



i{l + 2)q, 



e-i 



21 + 1 2f4 



2^ 2^ 



(2.59) 
(2.60) 
(2.61) 
(2.62) 



where = Za/r and is always less than unity for r outside the horizon. For the a solution, the multipole decomposition 



of the Wi vector, in the notation of Eq. (2.52), is 



(2.63) 



An 



21 + 1 4f 2 



(2.64) 



The general longitudinal solution for the extrins ic cu r vatur e, for each multipole, is a superposition of three contri- 
butions: (i) the particular solution given by Eqs. (p36|)-(^39|), (ii) the homo geneous "image" solution corresponding 



to Eqs. ( pl5| ), ( ^^ , for some ^-depend ent c hoic e of Zjmagc; multiplied by an arbitrary, ^-dependent constant, and (iii) 



the homogeneous "a" solution in Eqs. (2.57) - (2.64), for some ^-dependent choice of Za, multiplied by an arbitrary, 
^-dependent constant. One is led to ask what the "proper" choice is for the general longitudinal solution. A somewhat 
broader question is whether the initial data should be chosen to be longitudinal. 

In the initial data solutions generated for use in numerical relativity, the choice has often been made to take the 
solution to be inversion symmetric |lj]. (That solution turns out not, in fact, to be longitudinal.) The criterion 
used in numerical relativity has been definiteness and the simple-to-implement inner boundary condition provided 
by inversion symmetry. Here our criterion shall be quite different: we study what choice of hypersurface data is the 
best representation of the solution that evolves from an earlier astrophysical configuration. More specifically, we can 
start a particle falling from some large radius rg, and numerically evolve V' to a hypersurface at some later time. 
We can then study the agreement between that evolved data and the data given by some prescription for assigning 
hypersurface data to represent the falling particle. 

It should be understood that in making this comparison (i.e., evolved data vs. prescribed data) there are no issues 
of "choice" involving the slicing of spacetime. The Moncrief wave variable -0 is invariant with respect to perturbative 
changes in slicing (as well as shifts and diffeomorphism) , so that our comparison of evolved if) and prescribed ip involves 
no gauge related ambiguities. 



D. Choice of hypersurface data 

In our discussion above of initial data we found the following: (i) If we restrict ourselves to conformally flat data, we 
must choose a single constant, for each multipole, to specify our solution for ^; if we do not restrict to conformally flat 
data there is an unconstrained function of r for each initial multipole of ip. (ii) If we restrict our choice to longitudinal 
initial extrinsic curvature, then we must specify two constants for each multipole, in order to fix the initial ip; if we do 
not restrict to initial longitudinal extrinsic curvature, then there is an unconstrained function of r for each multipole 
of V' (in our case of axially symmetric collision) . 

Our goal below will be to see whether there is a simple physical choice that can be made for the prescribed 
hypersurface data that agrees well with the data evolved from t = Q. We shall try to reach that goal by restricting 
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our prescription to conformally-flat, longitudinal (hereafter CFL) data. The question we shall be asking, therefore, 
is whether there exists a choice of constants specifying the CFL solution for the particle on a hypersurface so that 
it agrees well with an evolved solution. Such a choice of constants could be determined by looking for some kind of 
"best fit" of the CFL solution to the evolved solution, but of much greater interest is a physically-based prescription 
for the choice of constants which can be used as a best guess of appropriate initial data when an evolved solution is 
not available. 

The rule we shall use for the specification of the CFL ip is based on a simple physical consideration. We fix the 
constant in the CFL solution for by demanding that the value of tp at the horizon be "frozen," i.e., fixed at the 
value it has on the initial hypersurface t = 0. This is in accord with the idea of dynamics freezing at the horizon, 
but it is better to view it as a consequence of the wave equation Eq. ( ^.4[ ), and its horizon limit at r* — > — cx). The 
initial data at t = is a smooth function of f near the horizon, and hence ip becomes constant as r* — )■ — oo. But the 
potential decreases as and can be ignored near the horizon, so ip = constant is a static solution of the wave 

equation. The value of ip at the horizon will therefore be unchanging in time. We shall see, in the numerical results 
presented below, that the evolved data demonstrate this freezing. 

It is straightforward to use this condition of freezing at the horizon to fix the constant multiplier of the homogeneous 



s oluti on in Eq. (2.2C), and the consequence is interesting. If we start with a BL solution (i.e., the solution in Eq. 
( ^.19 ) with the particle position at Zp = tq) then the "horizon frozen" solution, on a later hypersurface, turns out 
to be a superposition of the following simple solutions: (i) the BL solution corresponding to the the particle at its 
correct position Zp at time i, (ii) a solution for a particle inside the horizon located at ^imago — ('7i/2)^/zp, so that 
when added to the BL solution the total value of ip vanishes at the horizon, and (iii) a solution corresponding to a 
particle at the image location of the original position z = (M/2)^/fo. The magnitude of the mass parameter for this 
last image gives a ip at the horizon equal to the frozen value. 

The way in which we fix the constants for ^ is a more complicated, less elegant, implementation of the same basic 
physical idea of looking at the nature of solutions of Eq. ( p.4[ ), near the horizons. At i = our initial data is taken 
to be stationary, and of the form 

^ ^ a + bf + cf'^ -i . (2.65) 

Since f is related to r* by 

f = M/2 + A/e-l/2gr.V4M ^ Me-le'-*/4M ^ Q^^3r'/4M-^ ^ (2.66) 

it follows that our initial solution has the horizon limit 

■iP = h + fee"/-**^ + ge"/2*^ , ^ = . (2.67) 



In the horizon limit, the potential Vi, in Eq. (2.4), takes the form 



where Vo is a constant. In the horizon limit we are, then, solving the problem 



with the initial data given by Eq. ( |2.67 ). It is simple to show that the solution to this problem is 

ip=^a + 6cosh(</4Af)e'''/4*^ + [{d~Voi2Mfa) cosh{t/2M) + Voi2Mfa] e'^'/^M ^ (2.70) 

We can now take the time derivative of this near-horizon solution, evaluate it on the hypersurface at time t, and find 

^1, = {b/UI) sinh (t/4M)e'''/**^ + (S - Voi2Mfd.) sinh (t/2M)e''*/2*^ . (2.71) 

This result can be compared with the Taylor expansion around f ~ M/2 of the "basic," "image," and "a" longitudinal 
solutions, and the two constants fixing the solution thereby determined. 



From Eq. (2.70) it can be immediately seen that d does not depend on time. This is equivalent to t he co ndition of 



"freezing" on the horizon. The amplitude F by which we have to multiply the homogeneous solution (2.20) is 

i+i 

m2{zp) V^o 



'^homog[Zp) = ( ^ ) (1 -I- Tho,nog{ro)) - 1 , (2.72) 
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where Thomog{ro) represents the arbitrariness in choosing the data on the initial hypersurface t = 0. For instance, 
^homogiro) = 0, 1,-1 represents respectively BL, Misner, and antisymmetric data. 
The amplitude of the "image" term in the extrinsic curvature is the following 



r- — A 

^ image ^ 



M2 



6/(4; + 1) 

2/ + 3 



16mp(ro)(2/ + l)g|' 



4l{l~Thomogiro))sm{t/4M) 
4/3 _ 2/4^ 



(/ + l)(/2+/+l) 



■sin(i/2M) 



where 



M2 



2/ 



2/ + 3 



[-(5/ + l)(2/ + 3) + (/-6)(2/ + l 



Zimage = {M/2f/zp, Qs = M/{2zp), and 54 = M/(2fo). 
Finally, the amplitude of the a term in the extrinsic curvature (see Eq. ( ^.57j )) is 



-(/ + 7)(5/ + l) , (2/ + l)(26/2 + 26/-27) 2 

'I3 



(/-6)(5/ + 4) 



+ 



2/-1 (2/-l)(2/ + 3) 

16m2(ro)(l - ThomogiroWl + l)q'f\i^^ \ 5(/ - 2 



2/ - 1 



4m2(ro)(l + rhomogiro))q3 

/(/ + l)(/2+/+l)P 




10/^ + 9/ + 8 (/-6)(2/ + l) 2 
2/-1 + 



(2.73) 



(2.74) 



(2.75) 



where = (M/2)2/fo. 

Note that this choice of the amplitudes we assign to the homogeneous solutions fixes all three constants in the 
general solution, but this is not the only way we can fix them. We could have chosen, for instance, matching the 
behavior of the evolved for large r. This, however, leads to a less successful approximation. 



III. RESULTS 



A. Numerical method 



In this section we describe the algorithm used to integrate the wave equation (2.4) numerically. While the left hand 
side of the equation is straightforward to integrate, the source given by Eq. (2^) contains terms with a Dirac's delta 
and its derivative. Since we have not found in the literature a discussion of the numerical treatment of such sources, 
we shall describe the method in some detail. 

We have found it convenient to use a numerical scheme with step sizes Ai = ^Ar* = /i, and with a staggered grid. 
As Fig. H shows, this method connects points along lines of constant "retarded time" u = t — r* and "advanced time" 
V = t-\-r* . On this grid we have implemented a finite difference algorithm for evolving ■0 with errors of order h?' ; that 
is, at a given value of t and r* a measure of the solution error varies as h?' . This method has proven to be easy to 
implement and quite accurate. 



To derive our difference scheme we start by integrating Eq. (2.4) over the cell of our numerical grid shown in Fig. 
3. We use the notation 



dA 



dt dr* 



cell 



i+h 



du 



v+h 



dv 



(3.1) 



v — h 



Applied to the derivative terms in Eq. (p.4|) this gives: 



dA {~dtdt4> + ^ J J dA {-49„a„v} = 

4 [iP{t + h,r*)+ iP{t - h, r*) - V'(i, r* + h) - iP{t, r* - h)] 



(3.2) 
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Note that this resuh is exact; it contains no truncation errors. 

We next consider the integration of the potential term over the ceh. If the ceh is one with no source term contri- 
bution, then we can use 



dA {-Vip} = -h^ [V{r*)ij;{t + h,r*) + Vir*)'4j{t - h, r*) + 
V{r* + h)tp{t, r* +h) + V{r* - h)tP{t, r* ~ h)] + Oih"^) . 



(3.3) 



The /i^ order error in a generic cell is equivalent to an overall 0{h^) error in ip. 

The result in Eq. ( |3.3D assumes that ■tp is smooth in the grid cell. It cannot be applied to those cells through 
which the particle worldline passes, since ip is discontinuous across the worldline. For such cells we first obtain the 
coordinates (r*, ti) of the point where the particle enters the cell and (^2, t2) where the particle leaves it (see Fig. |^). 
Next, we divide the total area of the cell, (4/i^), into four subareas defined as follows: A2 is the part of the area of 
the rhomb below t — ti, A3 is the part of the area of the rhomb over t — t2, Ai is the remaining area to the left of 
the particle's trajectory, and A^ is the remaining area to the right. 

The integral of the Vijj term over the area of the cell is approximated by the sum of this function evaluated on the 
corners of the cell multiplied by the corresponding subarea Ai . This gives us 



dA {-V^} = -V{r*) [i;{t + h, r*)A3 + V(t, r* + h)A4 + ipit, r* - h)Ai + ^{t - h, r*)A2] 



(3.4) 



The truncation error in each such cell is of order (h^), just enough to have quadratic convergence, since only one cell 
with the particle has to be evaluated per time step. 

For cells through which the worldline passes, the integral of the source term in Eq. ( ^.8[ ) must be evaluated. As a 
convenience in discussing the numerical approximation of this term Eq. (^|^), we introduce the notation 

S = G{t, r)d[r - rp{t)] + F{t, r)d'[r - rp{t)] . (3.5) 

The integration over the cell, when done with due regard to the boundary terms generated by the S'[r — rp{t)], gives 

Ct2 



SdA = 2 



dt 



ti 



G{t,Tp{t)) 

1 - 2M/rp{t) 



— dr 



( F{t,r) \ 




V 1 - 2M/r J 





{l-2M/rp{h)f ^ ^ " 



±2- 



(3.6) 



The / dt integral in the first term can be performed to any precision since F and G are known functions. For our 
goal of quadratic convergence, a trapezoidal approximation for the integration is adequate. In the second term the 
upper (lower) sign is for particles entering the cell from the right (left), or equivalently for r\ > r* {rl < r*). In the 
same way, in the third term the upper (lower) sign is for particles leaving the cell to the right (left), or equivalently 
> r* (rj < r*). 

When the form of F and G, given in Eq. (|2.8D, are used in Eq. (pM the result is 



(1 - 2M/rpit)) 

G{t,rpit)) 
1 - 2M/rp{t) 



— dr 



Keo\rp-2M) 
(A + l)(Arp + 3M) ' 

Fit,r) 
1 - 2M/7 



eo = - 2Af/ro, 

Ke^'{l~2M/rp) 
' {X + l)rp{\rp + 3My 



[3M^ - iXMrp - A(A + l)r| - 6M(1 - el)rp] 



r-;{t) = -eo'^2M/rp~2M/ro 



(3.7) 



O ur n umerical s chem e, for cells through which the particle worldline does not pass, is to solve for iplt , r*), using 
Eq. (U) and Eq. (Q in the integral of Eq. (U). For cells containing the worldline, Eq. (jsj), Eq. (|3J), Eq. ( |3^ ) 
and Eq. (3^) are used. In summary, the evolution algorithm we use is 



iP{t + h, r*) = -V(i - h, r*) + [V'(i, r* + h)+ ^{t, r* - h)] 



-V{t*) 



(3.8) 
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for cells not crossed by the particle, and 



tpit + h, r*) = -tP{t-h,r*) 
+il;{t, r* -h) 



1 



4 

Vir*) 



{Ai + A3 



+ iP{t,r* +h) 



Vir*) 



(3.9) 



for the cells that the particle does cross. 

The above equations cannot, however, be used to initiate the evolution off the first hypersurface. If denotes the 
time at which we have the initial data, we lack the values V'(^o ~ h)^ necessary to apply the evolution algorithm. We 
can, however, use a Taylor expansion to write 



^ito -h,r*)^ ^{to + h,r*)- 2hdtHto, r*) + 0{h^) 



(3.10) 



The right hand side can be used in place of ^'(^o — h,r*) in the application of the algorithm to evolve off the first 
hypersurface. It is important to note that this substitution is valid only if 4'{t,r*) is not singular between t ^ ~ h 
and t = to — h. This requires that the particle worldline not cross the vertical line at r* between t = to — h and 
t = to + h. In setting up the computational grid, we have been careful always to avoid such a crossing. 

The numerical method used here, evolving initial data for a partial differential equation on a staggered grid, has 
little in common with the transform method used in Paper I, in which we studied only momentarily stationary initial 
conditions. A comparison of the two methods (in the case of momentarily stationary initial data) provided a powerful 
check of both methods as well as insights into the relative efficiency and accuracy of the methods. The agreement 
of the two approaches turned out to be excellent. For the goal of producing an evolved waveform, the numerical 
evolution method was found to be faster to a transform method by orders of magnitude, and to give more accurate 
results. 



B. Numerical results 



In this section we present computed results for the infall of a particle, starting from rest at initial coordinate 
position ro- The particle trajectory is a radial geodesic as described by Eq. (2.7). By "t — 0," we shall always mean 
the hypersurface at which the particle was at rest. We shall consider two values of ro/2Af : both 15 and 1.5. The former 
is a starting point where the influence of the background hole is reasonably weak. (The gravitational redshift, for 
example, is around only 3%.) This represents, then, an astrophysical starting point for which a Newtonian description 
would be a good approximation. The second choice, ro/2M = 1.5 is not a reasonable astrophysical starting point, 
but complements the first choice of ro, magnifies certain effects, and is useful for exploration and illustration. 

We will be interested in Cauchy data, both tp and tp, on subsequent t — constant hypersurfaces. (We repeat 
here an important feature of the formalism: the gauge invariant information tp and tp is unaffected by a change in 
slicing. When we say, therefore, a i = constant hypersurface, we mean only t = constant to zero order in the particle 
mass.) There are several ways in which we can specify which later slice we are considering, (i) We could specify 
the Schwarzschild coordinate time t for this hypersurface. (ii) We could specify the location of the infalling particle 
corresponding to time t; for this purpose we use the notation rp (or r*, the r* equivalent). Note that Vp is related to 
the position parameter Zp by rp = Zp{l + M/2zp)'^. (iii) We could specify the parameter P describing the particle's 
momentum at time t. (This is the P parameter that enters into the Bowen-York description of extrinsic curvature; 
see the discussion in Sec. II relating this parameter to the particle's 4-momentum components. A picture of a sliced 
Schwarzschild spacetime, illustrating the particle trajectory and the labeling of the spacetime, is shown in Fig. ^. 

In our results the crucial concept is the difference between "evolved" data and "prescribed" data. By "evolved" 
values of ip and -0 we shall mean the values that are found on some hypersurface after numerical evolution forward in 
time from the t = original hypersurface. By "prescribed data" we shall mean data that is chosen according to one 
of the prescriptions of Sec. II. Prescribed data, for our purposes here, is always chosen from the set of possibilities 
that we describe as conformally flat, longitudinal. This means that the 3-geometry is conformally flat, and the 
extrinsic curvature is "longitudinal," in the sense defined by York For ip, specific choices that we have described 
in Sec. 11 include: (i) BL data, the simplest choice for -0, (ii) Misner data, data that is symmetric with respect to 
reflection under f (M/2)^/f in the conformally related flat space, (iii) "frozen" data, the conformally flat solution 
of the hamiltonian constraint that has the value of ip at the horizon flxed to be the horizon value on the original 
{t = 0) hypersurface. Prescribed choices for ip include (i) the basic Bowen-York [Q solution, the simplest longitudinal 
solution, and (ii) "horizon matched" data, the longitudinal solution with the form of ^ (more specifically, the first 
and second derivatives with respect to r*) matched to the evolved solution at the horizon limit. 



14 



We first present comparisons of radiated energy. In Fig. |4| we show the quadrupole energy E2 radiated during the 
infall of the particle from tq = 1.5(2M). On the original t = hypersurface, the form of i/j is taken to be the Misner 
solution. (Since the particle is momentarily stationary at t — 0, we have ^/i = on this hypersurface, of course.) The 
plot shows the energy radiated from the particle for times after the particle is at a location r* = r*. These energies 
are computed from the solution that evolves from the data on the constant time hypersurface labeled r* . The "true" 
data are the data that evolved from the t = original hypersurface on which the particle was momentarily stationary. 
The energy computed in this way is the "true" total quadrupole energy emitted during infall from 15(2M), and has 
the value 8.1 x 10-^ml/{2M). 

The plot shows the energy generated when the true (i.e., evolved) and ^p, on a hypersurface, are replaced with 
prescribed data appropriate on that hypersurface for the position and momentum of the geodesically falling particle. 
To find the points for the dotted curve, for example, on each hypersurface, position and momentum of the infalling 
particle were calculated and used to generate Misner data for ip, and Bowen-York data for ip, on that hypersurface. 
These prescribed data were then numerically evolved forward in time, and the resulting radiation to a distant observer 
was computed. This radiated energy is presented as a function of the hypersurface label r* . Similarly, the solid curve 

shows the result of replacing the true hypersurface data with frozen data for ip and horizon-matched data for ip. The 
energies of the dashed curve are the result of replacing the ip data with horizon-matched data, and retaining the true 
ip data. 

These results show that the "frozen-matched" prescription are in slightly better agreement with the correct radiated 
energy than are the other choices. The choice of frozen-matched conditions is also more justifiable, since it is based 
on a physical consideration. Throughout the remainder of the paper, we shall consider the frozen-matched choice as 
the best in the set of conformally flat, longitudinal possibilities. 

In Fig. ^ what is more important than the comparison of prescriptions for hypersurface data is the fact that none 
of them is very good. Figure ^ gives some insight into the details of the failure. Here, for infall from rg — 1.5(2A/), 
the data for ip on each hypersurface, are replaced by frozen data; the true data for ^p (i.e., the data evolved from 
t = 0) are retained. It is important and interesting that the replacement of the tp data introduces a much smaller 
error than the replacement of ijj. This is a similar conclusion that one finds in the results of Baker et al. |l^ in a 
rather different context. 

The previous energy results, of course, are for infall from an astrophysically unreasonable initial radius. In Fig. 
^, energy results are shown for infall from tq — 15(2M), in which case the true total quadrupole radiated is 1.64 x 
10~^mo/(2M). Results are shown for three choices of prescribed hypersurface data. Again we see that the agreement 
with the true energy is reasonably good if the true (i.e., evolved) data for ip is retained. For the other cases the 
agreement is good only if the data is replaced on a hypersurface well before the particle reaches the peak of the 
potential Vf at around 1.55(2M). 

We ask next what the appearance of hypersurface data is, and how the true and prescribed forms differ. In Fig. |^ 
we show the form of the true quadrupole tp generated by a particle falling from rest at ro = 15(2Af ), with Misner data 
initially specified. The figure shows the manner in which ip{f*) evolves from its initial Misner form at i = 0. As the 
particle moves inward, the simple prescribed form of a single peak evolves into more complex shapes. Figure ^ shows 
that on late hypersurfaces, after the particle, a.i t ~ 100, r* w has passed through the region in which the Zerilli 
potential is strong, ip clearly contains the shape of in- and outgoing waves with the profile of quasinormal ringing. 

Figures ^ - |ll| compare the true ^ on a hypersurface with prescribed horizon frozen data on the same hypersurface. 
It should be noted that the magnitude of the discontinuity in the r derivative of "0 is related, through the hamiltonian 
constraint, to the mass of the particle. It must, therefore, be the same for the pairs of curves in these figures, and 
for all forms of -0 satisfying the hamiltonian constraint. Results are shown for several different hypersurfaces, and the 
meaning is clear. The prescribed ip never contains the complexity of shape that indicates the presence of radiation. 

In Fig. |l^, Ip is shown for Misner prescribed data on a sequence of hypersurfaces, this time for infall from ro = 
1.5(2Af). In Fig. ^ the frozen prescribed ip solution is shown on the same hypersurfaces. The comparison of the two 
figures shows the difference induced by the freezing of the value of ip at the horizon. 

The differences in the evolved and prescribed forms of ip suggest that at late times the form of ip is not conformally 
fiat initial data. This is explicitly demonstrated and quantified in Figs, ^and |l5|, in which the gauge invariant index 
of conformality, developed in Sec. II. B is plotted. The traveling bumps, with the appearance of radiation, confirm 
that the wave content of the true initial data is associated with its failure to be conformally fiat. 

We next look at results for ip. Figure |l^, |l^, |l^ compare the true (i.e., evolved) form of the quadrupolar ip with 
the horizon matched prescription, for infall from rest (and Misner data) at t = 0, ro/2M = 15. As would be expected, 
the difference between the evolved and the prescribed forms grows with time. 

The same comparison of true and horizon-matched ip is given in Fig. ^ and this time for infall from ro/2A/ = 
1.5(2Af). In all the comparisons of ip, the horizon matched choice was made for the prescribed data. The specific 



15 



choice, however, makes httle difference in the comparisons. In Fig. ^ a comparison of the true ip and two prescriptions 
is given for infall from ro/2M = 1.5. (Infall from the smaU radius magnifies differences of evolution and prescription.) 
The difference between the basic Bowen-York [|lj] prescription for tp, and the horizon matched prescription, is small, 
while the difference between either of them and the correct data for tp is very large. The horizon-matched prescription 
is a slightly better approximation than the "basic" one. 

Another way of looking at the difference between the true data and prescribed data is to investigate the effect the 
difference has on outgoing radiation. Figures Fig. 23, 2^ show the waveforms (that is, ip{t) at large r) generated 



by different types of prescribed data. In each case, the solid curve, labeled t — 0.0, shows the "true" waveform, 
the waveform generated by evolution of the original, momentarily stationary data at t = 0. This true waveform is 
compared with the waveforms evolved from prescribed data placed on later hypersurfaces. The figures show that 
prescribed data on the t/2M — 63.9 hypersurface leads to an outgoing waveform in reasonably good agreement with 
the true waveform. For prescribed data on t/2M — 93.6 (with the particle at r*/2M = 3.91, or rp/2M = 3.15), 
however, large differences are evident between waveform and the true waveform. The disagreement is less severe when 
the evolved data for ijj is retained. 

An analogous comparison of waveforms is shown in Fig. ^for infall from ro/2Af = 1.5. Here the effect of prescribed 
data on the phase of the wave is much more evident than in the ro/2M — 15 case. 

In the work of Abrahams and Cook p|, one of the motivations for the present work, the "close limit" was used to 
evolve prescribed data. This technique ||l2| , |l5| is applicable to a hypersurface at a time late enough that the colliding 
bodies can be considered to be inside a single nearly spherical horizon. In this technique one uses only the large-r 
parts of the initial data of the colliding bodies. This technique was found to be surprisingly successful in dealing with 
head on collisions of equal mass holes. We saw in Paper I, however, that for very unequal mass holes the close limit is 
valid only at extremely small separations. For comparison, we have applied the close limit to the present problem. On 
a sequence of hypersurfaces we replace the hypersurface data by the large-r form of the data, and evolve the results. 
The resulting energies are shown in Fig. |2^. These results, like those in Paper I, show that the close-limit calculations 
of energy are of limited usefulness. 



IV. DISCUSSION AND CONCLUSIONS 



The numerical results of the previous section give a range of comparisons between astrophysically evolved, and 
"prescribed" hypersurface data (for quadrupole modes.) The first lesson to be learned from these results is that for the 
particle limit the standard prescriptions for hypersurface data are not adequate for describing an astrophysical strong 
gravitational field. For a particle falling from rg — 15 (2M), for example, an approximation for the radiation using 
prescribed data on a late hypersurface seriously overestimates the radiated energy. Figure || shows that prescribed 
data for the particle at Vp = 3.6M (equivalent to r*/2M = 1.6) gives a radiated energy too large by a factor of 10. 
Even for a hypersurface with the particle in the marginally strong field at Vp = 5.7M (equivalent to r*/2M = 3.4) 
the energy is overestimated by a factor of two. 

This is in marked contrast to the results of Abrahams and Cook Q for the head on collision of equal mass holes. 
They prescribed standard data on a late time hypersurface and found predictions of energy in excellent agreement with 
predictions from numerical relativity. Baker and Li Q], looking further into this problem, found that the predicted 
energy was remarkably insensitive to the choice of hypersurface. Good predictions could be made with data specified 
on hypersurfaces over a fairly wide range of times. There are certain technical differences between the nature of 
the prescribed data we use here, and that used by Abrahams and Cook. In particular, they used fully symmetrized 
data, whereas our preferred choice of prescribed data is data that is conformally flat, longitudinal data that is frozen 
and matched at the horizon. This is surely not the origin of the different conclusions. For one thing switching 
from symmetrized hypersurface data to antisymmetrizcd, in numerical relativity computations, makes only a minor 
difference in the results. More important, the work of Baker and Li does not use symmetrized data, and the study 
by Baker et al. finds excellent results when symmetrized data is replaced by unsymmetrizcd approximate data. 
As a further check of this point, we have redone several of our computations using symmetrized prescribed data and, 
as expected, found only minor differences from the results for other prescriptions. 

The crucial difference between the failure of prescribed data here, and its success in the previous studies, must lie 
in the difference between the head on collision of equal mass throats, and of a large mass throat with a small mass 
("particle") one. It is not difficult to see why this should be. It is the same reason the "close limit" approximation 
is not successful in the collision of very unequal masses. (See Fig. ^ and Paper I.) For the collision of equal masses, 
a common horizon engulfs both colliding throats as they begin to get into each other's strong field influence. The 
generation of the large amplitude gravitational waves (i.e., the excitation of quasinormal ringing) more-or-less coincides 
with the disappearance of the individual throats inside a single connected horizon. For a hypersurface corresponding 
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to both throats inside a single horizon, it is plausible (and is found to be valid by the success of the calculations) 
that details of the geometry near the throats is not important since they cannot influence the region exterior to the 
horizon. For a collision of high and low mass throats the situation is very different. Strong radiation is generated when 
the particle is near a coordinate distance r = 1.5(2M), near the peak of the Zcrilli potential. For our t = constant 
slicing this corresponds to a hypersurface for which there is no common horizon, i.e., the particle is still well outside 
the horizon of the background hole. For this slicing, in fact, the particle always stays outside the background horizon, 
and the details of its local field are always in causal connection with infinity. The fault would not seem to be in the 
slicing. A t = constant slicing was the choice in the method used by Abrahams and Cook |^ and by Baker and Li [|j , 
and is an obviously natural slicing for the particle infall problem. 

We conclude from these considerations that the fault lies in the nature of the prescribed data we have been using, 
i.e., the standard choice of conformally flat longitudinal hypersurface data. As shown in the figures of the previous 
section, there is a dramatic difference between the true hypersurface data on a late hypersurface, and prescribed data. 
There are two strong reasons to suspect that the most important difference between true hypersurface data and any 
of the prescribed data is the difference in ^p. First, our argument above suggests that the failure of prescribed data 
in the particle limit lies in the failure to describe fields near the particle. In Figs. ^ - we sec that there is (at 
least visually) a large difference between the true and the prescribed data close to the particle. This difference can 
be ascribed to the constraints on conformally flat longitudinal data. The sharp variations in ?A(r*) are, intuitively, 
features that should couple strongly to radiation. The inability of the prescribed data to model these features, must 
be viewed as potentially important. The second reason for focusing attention on ip is that the hypersurface data for 
"0 appears to be much more important than that for ■0. This was the conclusion in the study by Baker et al. JTot . In 
that work, which used colliding throats, and no stress energy, Einstein's equations were linearized in the momentum 
of the throats and a clear identification could be made of how much of the radiated energy could be ascribed to the 
initial 3-geometry, how much to the momentum, i.e., the initial extrinsic curvature, and how much to the interaction 
of the two. It was found that except for extremely small initial momentum the radiation was almost completely due 
to the extrinsic curvature. This was called "momentum dominance" by Baker et al. Here we cannot make such a 
clear distinction. Due to the moving particle source, treated as a stress energy contribution, we cannot ascribe the 
radiation to the two kinds of initial data information. Nevertheless, the energy results in our Figs. § - ^ show that 
some other form of "momentum dominance" applies to the particle infall; the energy radiated is much more sensitive 
to the details of the extrinsic curvature than to the details of the 3-geometry. 

The insights provided by the particle limit have helped to clarify what is needed in the way of hypersurface data 
for numerical evolution. A direct consequence of the present study is the realization that prescribed data will not 
work as successfully for unequal mass head on collisions as for the equal mass case. At some ratio of masses of the 
colliding holes the use of prescribed data on a late time hypersurface will start to give a significant overestimate of 
the radiated energy. We have, so far, studied only the quadrupole perturbations. Higher multipole moments are less 
important astrophysically, and are not likely to lead to very different conclusions. 

The head on collision, of course, is neither astrophysically plausible, nor an interesting source of outgoing radiation. 
The interesting case is the last stage in the orbital decay of a binary pair of holes. For this one would like to start 
with astrophysically reasonable data on a "as late as possible" hypersurface. An understanding of how to do this for 
the head on case is a necessary step (and perhaps a sufficient one) towards an understanding of the more general 
problem. With that motivation, we shall, in a subsequent paper, investigate what can be done to provide good late 
time data. We shall, in particular, abandon the traditional choice of conformally fiat, longitudinal data, and, instead, 
shall focus on data that gives a good description of the field near the particle. 
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FIG. 1. Geometry in the conformally related flat space. 

FIG. 2. A cell of the computational grid containing a segment of the particle worldline. Grid nodes are shown coordinatized 
with t, r* labels, and with values of the null coordinates u = t — r* and v = t + r*. The areas Ai - ■ ■ A4 are used as weights in 
the numerical algorithm. 

FIG. 3. Hypersurfaces of constant t in the Schwarzschild spacetime. The trajectory is shown of a particle starting from 

rest at position ro/2M — 15, and falling inward on a radial geodesic. Also shown arc hypersurfaces of constant time at 
t/2M = 20,40,60 and 80. The hypersurfaces arc labeled with the value of t, of the particle position Vp, and of momentum 
parameter P. The t/2M = 20 hypcrsurfacc, for example, can also be referred to as the P/mo = 0.05 hypersurface, the 
rp/2M = 14.6 hypersurface, or the r*/2M = 17.2 hypersurface. 

FIG. 4. Radiated Energy, for infall from ro = 1.5(2M) Energy is shown as a function of the r* label of a hypersurface and 
results are given for several types of prescribed hypersurface data. Radiated energy computed is always the energy radiated 
starting from the time of the hypersurface on which data is prescribed. See text for details. 

FIG. 5. Radiated Energy, for infall from ro = 1.5(2M). Energy is shown as a function of the label of a hypersurface. On 
that hypersurface the true, i.e., evolved, data for tji has been retained, but frozen prescribed data for has been substituted 
for the true V'- 

FIG. 6. Radiated Energy, for infall from ro = 15(2Af ). Energy is shown as a function of the label of a hypersurface and 
results are given for several types of prescribed hypersurface data. The computed energy is always the energy radiated starting 
from the time of the hypersurface on which data is prescribed. See text for details. 

FIG. 7. The true quadrupole i/'(r*) for infall from ro = 15(2M), i.e., tp computed by numerical evolution from the initial 
hypersurface. Curves are shown for several hypersurfaces labeled with r^, the value of the particle's r* location on each 
hypersurface. 

FIG. 8. The true quadrupole ^p(r*) for infall from ro — 15(2Af), shown on late hypersurfaces. For t/2M = 109.7 the shape 
of tp shows outgoing quasinormal radiation; for t/2M = 149.6 both ingoing and outgoing quasinormal oscillations are evident. 

FIG. 9. The comparison of true and prescribed quadrupole ^ for infall from ro = 15(2M). The form of the true data tp{r*) 
for infall from Misner data on the original hypersurface, is compared with the frozen prescribed data on the hypersurface at 
t/2M = 52.2. 



FIG. 10. For infall from ro = 15(2M), the evolved data is compared with the frozen prescribed data, at t/2M = 96.9 
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FIG. 11. For infall from ro = 15(2M), the evolved data is compared with the frozen prescribed data, at t/2M = 100.2. 



FIG. 12. Misner prescribed data for infall from ro = 1.5(2M), shown on a sequence of hypersurfaces. 



FIG. 13. Frozen prescribed data for infall from ro = 1.5(2M), shown on a sequence of hypersurfaces. 

FIG. 14. The index of conformality /conf for infall from ro = 15(2M) with Misner initial data. The index, for quadrupole 
perturbations, is equivalent to H2—K in the RW gauge, and is a gauge invariant measure of the extent to which the hypersurface 

3-gcomctry fails to bo conformally flat (see text). The index is given, as a function of r* , for several difi'oront hypersurfaces. 
The noisy nature of the curves is caused by the need to take second differences of numerical results to compute /conf. 



FIG. 15. The index of conformality Jconf for infall from ro = 15(2M), shown on late hypersurfaces. 



FIG. 16. For ro/2M = 15, the true (i.e., evolved) form of tp compared with the horizon matched prescribed form of 'ip, for 
t/2M = 52.2. 



FIG. 17. For ro/2M = 15, the true (i.e., evolved) form of ip compared with the horizon matched prescribed form of ip, for 
t/2M = 96.9. 



FIG. 18. For ro/2M = 15, the true (i.e., evolved) form of ip compared with the horizon matched prescribed form of ip, for 
t/2M = 100.2. 



FIG. 19. For ro/2M = 1.5, the true (i.e., evolved) form of ip compared with the horizon matched prescribed form of ip, for 
rj,/2M = 1.34. 



FIG. 20. For ro/2M = 1.5, the true (i.e., evolved) form of ip compared with the horizon matched prescribed form of ip, for 
rp/2M = 1.21. 

FIG. 21. For ro/2M = 1.5, the true, i.e., evolved, form of ip compared with the horizon matched prescribed form of ip, and 
with the Bowen-York form. 

FIG. 22. The quadrupole waveform, ip{t) at r/2M = 1000, for infall from ro/2M = 15. Three curves are shown corresponding 
to the true waveform (data prescribed at t/2M = 0) and to the waveforms generated when frozen ip data and horizon matched 
Ip data are imposed on later hypersurfaces. 

FIG. 23. The quadrupole waveform, ip{t) at r/2M = 1000, for infall from ro/2M = 15. Three curves are shown corresponding 
to the true waveform (data prescribed at t/2M = 0) and to the waveforms generated when horizon matched ip data are imposed 
on later hypersurfaces and the true, i.e., evolved, values of ip are retained. 

FIG. 24. The quadrupole waveform, ip{t) at r/2M = 1000, for infall from ro/2M = 15. Three curves are shown corresponding 
to the true waveform (prescribed data a t/2M = 0) and to the waveforms generated when frozen ip data are imposed on later 
hypersurfaces and the true, i.e., evolved, values of ip are retained. 
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FIG. 25. The quadrupole waveform, at r/2M = 1000, for infall from ro/2M = 1.5. The sohd curve shows the "true" 
waveform evolved directly from the initial data at t/2M = 0. The dashed curve gives the waveform generated when prescribed 
data (horizon-frozen data) for ip, and evolved data for ip are placed on a late hypersurface. 

FIG. 26. Radiated energy for infaU from ro = 15(2M) computed with the close limit. On every hypersurface close limit data 
was substituted and evolved. The computed energy is plotted as a function of the hypersurface parameter on which the close 
limit data was substituted. 



20 



ro=1-5(2M) 

0.8 I ^ 1 ' 1 ' r 




r*p /(2M) 



LU 

CM 

C 

E 

CM 



0.015 



ro=1-5(2M) 



0.010 



0.005 



frozen-evolved 



0.000 
-2.5 



-1.5 



-0.5 



0.5 



r* /(2M) 



ro=15(2M) 




frozen-matched 

evolved-matched 
frozen-evolved 



0.00 



9 

r* /(2M) 



13 



17 



ro=15(2M) 



2.8 



2.4 



2.0 



1.6 



1.2 



0.8 



0.4 



0.0 



-0.4 



-10 



10 







1 










\ 


r*p=17.6 , t=0.0 




- 


1 ii 

! ii 




r*p=7.2 , t=85.7 
r*p=2.1 , t=96.8 


- 








r*p=0.0 , t=99.5 




- 


1 ' ' 

' : 






- 


1 

1 

1 

/ 


' I 
' 1 

1 ' 1 
1 ' ' 

\ I J 










\ 









30 



50 
rV(2M) 



70 



90 



IK 



2.0 



ro=15(2M) 



1.5 



r*p=-10,t=109.7 
r*p=-50,t=149.6 



1.0 



0.5 



0.0 




-0.5 



-80 -40 



40 80 

rV(2M) 



120 160 



ro=15(2M), r =12.0(2M), P=0.14(mo), t=52.2(2M) 

3.0 




-20.0 20.0 60.0 100.0 

rV(2M) 



ro=15(2M), r =2.0(2M), P=1.27(mo), t=96.9(2M) 



frozen 
evolved 




60.0 

rV(2M) 



100.0 



ro=15(2M), rp=1.16(2M), P=5.0(mo), t=100.2(2M) 



3.2 



2.6 



2.0 



1.4 



0.8 



frozen 
evolved 




60.0 

rV(2M) 



100.0 



ro=1.5(2M), Misner 




-10 -5 5 

rV(2M) 



ro=1.5(2M), Frozen 



r^=1.50 




-10 



-5 







rV(2M) 



ro=15(2M), 1=2 



0.03 



r*p=17.2, t=20 
r* =15.8, t=40 



r* =13.2, t=60 



0.01 



-0.01 



-0.03 




-0.05 



-15.0 



15.0 



45.0 

r7(2M) 



75.0 



105.0 



0.20 



0.15 



0.10 



0.05 



0.00 



-0.05 



-0.10 



-20 



ro=15(2M), 1=2 



r*p=7.5, t=85 
r*p=5.6, t=90 
r*p=3.2, t=95 
r*p=-0.5, t=100 








20 



40 60 

rV(2M) 



80 100 120 



ro=15(2M), 



rp=12.0(2M), P=0.14(mo), t=52.2(2M) 



0.010 



-0.010 



o 

E 

CM 



-0.030 



-0.050 




evolved 
matched 



-0.070 



-25 







25 50 

rV(2M) 



75 



100 



0.5 



0.0 



^ -0.5 



o 

E 

CM 



-1.0 



-1.5 



-2.0 



-5 



ro=15(2M), r =2.0(2M), P=1.27(mo), t=96.9(2M) 







5 10 

rV(2M) 













^ / 

X / 

- 

- 


/ / 

/ ^ 

/ / 
/ / 
/ / 
/ / 




/ / 
/ / 
/ / 
/ / 


evolved 

matched 






/ / 
/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 



15 



20 



1.0 



ro=15(2M), rp=1.16(2M), P=5.0(mo), t=100.2(2M) 



-1.0 



CM 



-3.0 



evolved 
matched 



-5.0 



-10 







rV(2M) 



10 



20 



ro=1.5(2M), rp=1 .34(2M) 



0.2 



CM 



-0.2 



-0.8 




evolved 
matched 



-1.2 



-10 







r7(2M) 



10 



20 



ro=1.5(2M), r =1.06(2M) 



0.0 



o 

E 

CM 



-1.9 



-3.8 



-10 







/ 




^ \ 

- 


/ 

1 / 

./ 

/ 

;/ 
;7 

./ 




- 
- 


- 


:/ 

1 


evolved 

matched 
basic (BY) 






\ il 
\\ // 
\\ // 
\\ // 








v\ // 

V-7 







10 



20 



rV(2M) 



ro=15(2M), r^b3=1000(2M), 1=2, (frozen-matched) 



950 



r*p=17.6 


t= 


=0.0 , 


P= 


0.0 


r*p=12.6 


t= 


=63.9, 


P: 


=0.19 


r>3.91 


t= 


=93.6, 


R: 


=0.73 




1000 



1050 
t/(2M) 



1100 



1150 



ro=15(2M), r„b3=1000(2M), 1=2, (evolved-matched) 



0.50 



r*p=17.6 


t= 


=0.0 , 


P= 


0.0 


r*p=12.6 


t= 


=63.9, 


P: 


=0.19 


r>3.91 


t= 


=93.6, 


R: 


=0.73 



0.30 



0.10 



-0.10 



-0.30 
950.0 




1000.0 



1050.0 

t/(2M) 



1100.0 



1150 



ro=15(2M), r , =1000(2M), 1=2, (frozen-evolved) 



0.50 



r*p=17.6 


t= 


=0.0 , 


P= 


0.0 


r*p=12.6 


t= 


=63.9, 


P: 


=0.19 


r>3.91 


t= 


=93.6, 


R: 


=0.73 



0.30 



0.10 



-0.10 



-0.30 




950 



1000 



1050 
t/(2M) 



1100 



1150 



